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Abstract 

We propose a general physical mechanism which could contribute to the 
formation of fast line-driven outflows at the vicinity of strong gravitational 
field sources. We argue that the gradient of the gravitational potential plays 
the same role as the velocity gradient plays in Sobolev approximation. Both 
Doppler effect and gravitational redshifting are taken into account in Sobolev 
approximation. The radiation force becomes a function of the local velocity 
gradient and the gradient of the gravitational potential. The derived equa- 
tion of motion has a critical point that is different from that of CAK. A 
solution that is continuous through the singular point obtained numerically. 
A comparison with CAK theory is presented. It is shown that the developed 
theory predicts terminal velocities which are greater than those obtained 
from the CAK theory. 
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Introduction 

Acceleration of matter due to radiation pressure in lines plays an important 
role in formation of outflows from hot stars and possibly from active galactic 
nuclei (QSO, Seyfert galaxes, etc). This paper is the first in a series of papers 
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on the effects that strong gravitational fields play in formation and structure 
of winds driven by radiation pressure in lines. 

Theory of winds from O-type stars is well-developed and has a good 
agreement with observations. In a pioneering paper of Sobolev (1960), it 
was recognized that a radiation transfer in accelerated medium is simplified 
drastically in comparison with that of static atmosphere. The importance of 
the line opacity for the formation of winds from hot stars was pointed out in 
paper of Lucy & Solomon (1970). A prominent step in this field was made in 
papers of Castor, Abbott & Klein (1975), ( hereafter CAK ). In these papers 
it was shown that absorption of the radiation flux in lines can be an effective 
mechanism of redistributing of momentum from radiation to the wind. The 
presence of the velocity gradient gives additional effect on the acceleration 
because of the Sobolev effect. CAK showed that the resultant radiation force, 
which is due to absorption in an ensemble of optically thin and optically thick 
lines, may be several orders of magnitude greater than that due to electron 
scattering. The ideas and technique of CAK were developed by Abbott 
(1980) and in works of many other authors. 

Theory of winds accelerated by the radiation pressure in lines is usually 
applied to explain outflows from active galactic nuclei (AGN). The many 
puzzling, observational characteristic features of AGNs are: the blue-shifted, 
relative to the emission line rest frequency broad absorption lines (BALs) 
- the most convincing evidence of outflows with velocities as large as 0.1c; 
NALs -narrow absorption lines seen in UV and X-ray spectra associated 
with ~ 1000 km • s _1 outflows; BELs - broad emission lines observed in UV 
indicate the flow with the characteristic velocity about several thousand 
km-s" 1 . From the short time-scale of the X-Ray variability (Tennant et 
al 1981) it is concluded that the size of the emitting region is ~ 10 14 cm. 
(Although there exist an uncertainty in the estimation of the column densi- 
ties, see Arav et al, 2002, and references therein) Only black hole (BH) as 
a central object can couple this small length scale with the total luminos- 
ity L ~ 10 46 egr • s -1 . The existence of BALs together with large luminosity 
makes radiation driven mechanism plausible. Line-locking, observed in spec- 
tra of some QSO gives additional evidence of the importance of radiation in 
acceleration of the wind. Murray at al. (1995) model the wind that origi- 
nates not far from BH ~ 10 16 cm. Simplifying assumptions allowed to solve 
separately equation of motion in radial and polar angle directions. Two- 
dimensional time-dependent hydrodynamic calculations of radiation-driven 
winds from luminous accretion disks were made in Proga et al. (1998), Proga 
et al. (2000). 

In all papers where the radiation pressure was assumed as a main driv- 
ing force forming outflows from AGNs, authors adopted modifications of the 
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CAK theory (although effects which are not included in the CAK theory 
were taken into account). A serious difference between a hot star wind and 
AGN outflow is due to different geometry. Wind from O-type star is nearly 
spherically-symmetric whereas an outflow from AGN originates from accre- 
tion disc an thus is closely axially-symmetric. From observations we know 
that wind from AGN is exposed to a very high X-ray and UV flux. The AGN 
spectra is harder then that of O-type star thus the radiation force could be 
smaller due to a higher ionization state. On the other hand the observations 
of spectral lines from moderately and highly ionized species pose a problem 
of how the wind avoids over-ionization. This problem can be overcame if 
there exist some gas that shields the wind from the ionizing radiation. The 
second way is to assume that the wind is in the form of dense clouds which 
are confined by some mechanism. The ionizing effects of the X-ray flux in 
the winds from massive X-ray binary systems were considered in Stevens & 
Kallman (1990). 

The main underpinning assumption of CAK, is that of the velocity gra- 
dient leads to enhancement of the radiation force. If there is no velocity 
gradient then the radiation force is simply due to absorption of the fraction 
of the radiation flux which is blocked by lines. If there is a velocity gradient, 
the absorbing ions are exposed to the unattenuated radiation flux due to 
Doppler shifting. 

In this paper we propose a general physical mechanism which could play 
an important role in the formation of fast outflows at the vicinity of strong 
gravitational field sources. We argue that the gradient of the gravitational 
potential d(f>/dr (in case of sufficiently strong gravitational fields) plays the 
same role as dv/dr plays in acceleration of line-driven winds from hot stars. 

Taking into account gravitational redshift of frequency we conclude that 
the gradient of the gravitational potential allows a line to shift out from be- 
hind of its own shadow. The line is shifted to the extent where the continuum 
flux is not reduced by absorption, hence d<fi/dr exposes the line to powerful 
radiation. Thus we call the resultant wind: " Gravitationally Exposed Flow" 
(GEF). 

In this paper the following toy model is adopted. Consider a wind accel- 
erated by the radiation pressure in lines at the vicinity of BH. For simplicity 
we assume the wind to be spherically-symmetric and skip ionization prob- 
lem. These crude assumptions are justified by the question we are going to 
answer that is to compare a solution for gravitationally exposed flow with 
CAK solution. We are addressing 2d modeling of the problem to the separate 
paper. 

The plan of the paper is as follows. In Section 1. we discuss relations 
for the optical depth r L and radiation force f r in standard CAK theory. 
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In Section 2. a relation for the radiation force is derived taking into ac- 
count gravitational redshifting. Equation of motion describing stationary, 
spherically-symmetric wind is derived. Properties of the solution at the crit- 
ical point are investigated. The numerical solution is provided. The results 
are summarized in the Discussion and Conclusions. 



1 Radiation force 

Assume that a radiation flux F c (u) (F c (ui) - is the continuum flux per interval 
of frequency on the line frequency Ui) emitted by the photosphere (or disc) 
is radially-directed. The radiation force that results from absorption in a 
single line with the Doppler width Avd ~ vvth/c may be represented in the 
following form (Castor, Abbott & Klein 1975): 

fl = Kl (^fM) IzfL, (1) 

The effective optical depth t\ is calculated on the frequency of line, which 
obtained taking into account the Doppler shifting, k\ (cm 2 g _1 ) is the opacity 
at the line centre: 

Tie 2 gf Ni/gi - N u /g u 

K i * • \ z ) 

mc p Ave) 

In Sobolev approach dv/dr is taken into account when calculating 77. A 
relation for the optical depth was derived by Castor (1970): 

ti= f ippnidr ~ r! th !\ L , ■ (3) 
Jr \dv/dr\ 

When calculating 77 it is necessary to take into account only those absorbers 
which are in a section of the column across which the velocity changes by 
v t h- Relation @ was obtained assuming that the line profile ip is a delta 
function. CAK introduced an optical depth parameter independent of line 
opacity: 

t = Vep-p^, (4) 
dv/dr 

where a e = 0.4 cm 2 ■ g -1 . If the velocity gradient is large the optical depth 
parameter t may be much less then the corresponding electron optical depth. 
If t ~ 1 the radiation force is determined by the part of the radiation flux that 
is blocked by lines. In the opposite case of small t, the radiation force may 
exceed gravity by several orders of magnitude for a typical O-star atmosphere. 
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Summing equation ([]]) over an ensemble of lines, CAK obtained the following 
relation: 

!> = ^-Mif), (5) 
where M(t) is the force multiplier: 

M{t) = Y. FM f VD ™<^t-\ (6) 

i 

where (3l = ^l/^c (note that t = tl/(3l), F - is the total radiation flux. 
Ti/ 'ki - is the column of mass which absorbs radiation, (1 — e n ) multiplier in 
([!]) gives the probability of absorption. If there are only strong lines then 
the radiation force is proportional to the number of strong lines. Thus in 
optically thick case f\ is independent of the line strength and fi ~ dv/dr. 
If lines are all optically thin then each ion is absorbing the unattenuated 
radiation flux, and thus the radiative acceleration is independent of the wind 
dynamics. For the ensemble of optically thin and thick lines CAK found that 
M(t) can be fitted by the power-law: 

M(t) = kt~ a . (7) 

a = 1 represents when all lines are optically thick, a = - optically 

thin. With (|?|) and (f|), equation (|5|) reads: 

La e k ^ a (^dv\ a 

h = - — -{(TeVth) \--r) ■ (8) 
Airr z c \p dr J 



2 Gravitationally Exposed Flow 

2.1 The radiation force due to gradient of gravitational 
field 

If a wind is accelerated in a strong gravitational field we may expect that a 
photon emitted by the disc may become resonant with the opacity of some 
line not only because of the Doppler effect but also because of the grav- 
itational redshifting. The redshifting of the photon's frequency is the only 
effect of general relativity (hereafter GR) we take into account in our approx- 
imate treatment. A photon with the frequency z/^, emitted at the point with 
the gravitational potential <p\ will be registered at the point where = 02 
(| 0i |>| 02 I, < 0) reddened according to the well-known approximate 
formulae (see for example Landau Sz Lifshitz, 1960): 
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If there is a velocity difference between these two points the photon will be 
additionally red-shifted due to the Doppler effect. The resultant frequency, 
seen by the absorber at his rest frame reads: 

v*u i+ u d (-z+wy (io) 

To derive a relation for the radiation force we need to calculate the cor- 
responding Sobolev optical depth, but now taking into account the grav- 
itational redshifting. Assuming that absorption occurs in a line with the 
delta-function profile, for the optical depth in a radial direction we obtain: 

f°° , PVth^i nr , 

Tl = l KlpdT ^ dv/dr+^/dr - 

A characteristic length: 5r ~ v t h/(dv/dr + -d(p/dr) gives a thickness of a 
shell where the absorption due to a single line take place. 

In our calculations we are neglecting the influence of the gravitational 
field on the geometry of space-time, because in that case it our treatment 
will be too complicated. Introducing an optical depth parameter, which is 
analogous to (f|), we have: 

t _ °epvth , 12 * 
~~ dv/dr + \d(j)/dr' 

Note that relations @, (0) are derived in such a way that all information 
about wind dynamics (distribution of v that is responsible of redshifting) is 
imbedded into parameter t. It is assumed that r must be calculated tak- 
ing into account redshifting. No assumptions about the particular physical 
mechanism that produces redshifting are made. In our case t contains also an 
information about additional redshifting which is due to gradient of gravita- 
tional potential. Note that (|12D may be represented in the form t ' 1 



dv*/ dr ' 

where v* — v H — <fi. The distribution of v* contains all the information about 

c 

redshifting. These considerations hold so far Sobolev approximation is as- 
sumed and until the influence of the gravitational field on the geometry of 



space-time is neglected. We are justified in making use of fll2[) instead of 
OD) when calculating (^). 



By substituting (fl2f) into © and (0) the radiation force is obtained: 
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vr 



' dv ld(f) 
dr c dr 



(13) 



]M\av th/ 

where the radiation force was transformed making use of the continuity equa- 
tion: 



M |= Aixpvr 2 



(14) 



2.2 Basic equations 

An equation of motion describing stationary, spherically-symmetric, isother- 
mal, line-driven wind reads: 



dv 

! 

dr 



1 dP dip La La e k 



47T 



p dr dr Airr 2 c Airr 2 c \ \M\av t h 



vr 




(15) 



As well as the equation of motion of CAK, equation (|1|) is nonlinear with 
respect to the velocity gradient. This bahaviour complicates a structure 
of (|15D and slows numerical solution. Investigating (|15D we hope to find a 
trans-critical solution, which starts subsonically, proceeds through the critical 
point and goes to infinity approaching a terminal velocity v°°. Equation 



( lq) resembles CAK equation of motion except the radiation pressure term. 
Instead of (dv/dr) a term, equation (|15|) includes a combination: (dv/dr + 
-d(p/dr) a . As we will see this property of the equation (|T5| ) changes the 
structure of transcritical solution. 

We adopt two types of potentials: Newtonian potential (NP) and Paczyn- 
ski and Wiita (hereafter PW ) potential: 



GM 



(16) 



We make use of the modified potentials to model approximately the GR ef- 
fects. PW potential, for the Schwarzschild black hole correctly reproduces the 
positions of both the last stable orbit and the marginally bound orbit. The 
many properties of modified potentials are discussed in Artemova, Bjornsson 
& Novikov (1996). Introduction of the modified potentials allows to take 
approximately into account some properties of the exact relativistic descrip- 
tion. For simplicity hereafter we use Newtonian potential. The case of the 
PW-potential is described in Appendix. For Newtonian gravity <fi = —GM/r, 
equation (|15D reads: 
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4irr 2 c \ \M\av th 
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dv 
dr 



4 7rr 2 c 
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c r' 



where a = \ —— ] 

was assumed that the wind is isothermal T 
Introducing nondimmensional variables: 



= 0, (17) 

a sound velocity. When obtaining equation ( ]T7| ) it 



const, and P = a 2 p. 



r 

x = — , 



v = — , 

V c 



(18) 

where r c is the critical point radius and v c is the velocity of matter at the 
critical point. Taking into account (|T8|), equation (|TTD reads: 



F(p, v, x) 



vx 2 p 



2a\x + y (1 







r) 



(19) 



where p 



— and T 

dx 



For simplicity hereafter we omit tilde. The 
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following nondimmentional parameters were introduced: 



at 



a 



c 

V r . 



(20) 



Note that j3 = —a\ and there is only one independent parameter a\. A con- 
a 

stant /i is determined according to relation: 



H = Tk 



8ttGM 
M I a e v th 



(21) 



Equation ( Jl9| ) is non linear with respect to p. The point where the speed 
of sound is equal to velocity of the flow is no longer the singular point of 
the equation of motion. To treat such equations a special technique must 
be applied. We are interested of the solution that starts subsonic near BH 
is continuous through the singular point and goes supersonically to infinity. 
Singular point is defined by the condition: 



8 



The second condition comes from the fact that at the critical point the ve- 
locity gradient is continuous which requires that v" is defined at the critical 
point. The regularity condition reads: 



'dF\ 



OF s 



-5T +Pc\-Z- 



dx J 



dv 



0. 



(23) 



where p c = — ( — j . 

v c \drj c 

From fll^D, ( p2|) the following useful relation is obtained: 

" 9 / n \ 1 1—Ot 

xi 



fi = a x (l — a\) 



Pc + 



2x1 



where x c = — . Substituting ( p!9|) to (p2|) , taking into account 
result in the following equation: 



(24) 
will 



2a{xy c + -{d{ - l)Pp c - 2a{ = 0. 



(25) 



The momentum equation (|l^) is valid a hole domain of the investigation, 
thus calculating it at the critical point and substituting /i from (|24"D will 
obtain the following equation: 



1 (1 - a\)x 2 c p c - 2a\x c + ? [ ' 



T) 



2^ P 



0. 



(26) 



a 2 17 2a 

Equations (p5|) , (|26|) form the system of equations from which p c and a x may 
be calculated. For example, assuming that the critical point is situated at 
100 r g , from (f25|), @) we find: a x = 0.0393124, p c = 10.4831. As it will be 
clear from the following, equations (|25|) - (|26|) are not only more complicated 
than those of CAK, but also manifest the different physical behavior of GEF. 
Equation (^) differs from those of CAK by the second term. Dropping 
in (E5|) this term will reveal the result of CAK for the isothermal wind: 



p c = ±l/xl, or f ^ 
into account that /3 - 



±— . To solve 



for ai, Pc, we must take 



Vth 



-ai, where v t h 



T 
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2.3 Numerical solution 



For a given position of the critical point x c = r c /r g equations fl2"5|), ( |26"D are 
used to obtain the velocity gradient at the critical point. In contrast to CAK 
wind, the GEF case is more difficult to analyze. In the isothermal limit of the 
CAK wind, equations (p5|) , (f26|) could be solved analytically. Anfortunately 
it is not possible (albeit see the end of this section) in GEF case. In order 
to obtain p c and a± equations (p5|), (|26|) must be solved numerically. All 
calculations in this paper are made for 1O 7 M black hole, the temperature 
of the flow is T = 4 • 10 4 K other parameters are: V = 0.5, a = 0.6, k = 0.1. 

In case of stellar winds a photospheric boundary condition is usually 
adopted. Adjusting the position of the critical point one obtains a solution 
that gives the position of the photosphere (r = r p h at r e = 2/3) at some 
prescribed radius R, which is identified with the radius of a star (Bisnovatyi- 
Kogan 2001). Similar procedure was adopted by CAK. In case of AGN 
such "photospheric" conditions are clearly unphysical. The solution for the 
wind should be continuously fitted with the solution for the accretion disc. 
This requires 2D modelling which is beyond the scope of this paper. In 
spherically-symmetric approximation adopted here it is not possible to fit 
self-consistently a wind solution with that of accretion disc. 

To compare GEF solution with that of CAK we start deeply subsonic 
(v <C v s ), starting from some initial density pi n and calculating both CAK 
and GEF solutions. Mathematically it is equivalent to the problem of fitting 
of the wind solution with that of a static core when calculating a structure 
of a star when mass loss is taken into account. In such solution 
for the outflowing envelope is continuously fitted with that of a static core. 
As it results from a stellar wind theory (for stationary spherically-symmetric 
wind) to fit continuously a solution for a stationary outflowing envelope with 
that of a static core, only p and T should be fitted. In case of stellar wind 
this condition reads: p e ™ = p c °^t ■> = T%% e where the fitting point must 
lay in a deep subsonic region t)Ci; s (see Bisnovatyi-Kogan & Dorodnitsyn 
1999 and references therein). 

The adopted procedure is equivalent to that if we start from some static 
configuration (v << v s ) with the density pi n and then relax to stationary 
wind solution (CAK and GEF). Prescribing p in at the very bottom of the 
wind (v <C v s ) at r in we may find the solution which satisfies boundary 
condition. This requirement allows to compare GEF solution with that of 
obtained from CAK theory. 

We step out from the critical point x c by means of the approximate formulas: 
v = 1± | p c | Ax, where Ax = x — x c . On every integration step equation 
(|j~9| ) is solved numerically in order to obtain dv/dx. The forth order, adaptive 
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Runge-Kutta scheme was used to obtain solutions depicted on Figure 1. - 
Figure 4. 

Adjusting x^ AK , x^ EF we select those solutions which both satisfy in- 
ner density condition. The following densities at the inner boundaries were 
adopted: p in = 10~ 12 gcm~ 3 for x c = 50, p in = 10~ 14 gcm -3 for x c = 100, 
Pi n = 10~ 16 gem -3 for x c = 500, pi n = 10~ 21 gcm~ 3 for x c = 8000. The com- 
parative results of the numerical integration of (|19|) for s = 1 (GEF solution) 
and for s = (CAK wind theory) are shown on Fig.l. These solutions were 
obtained for the following set of parameters: 



„GEF 


„CAK 


oo cm 
U GEF s 


„,oo cm 
V CAK T~ 


M GEF f 


M CAK f 


50 


75 


5- 10 9 


3.67- 10 9 


2.56597- 10 26 


2.55692 • 10 26 


100 


155.2 


3.3 • 10 9 


2.6 • 10 9 


2.56597- 10 26 


2.55693 • 10 26 


500 


751.7 


1.37- 10 9 


1.15 • 10 9 


2.56599 • 10 26 


2.55697- 10 26 



The difference between GEF and CAK wind is more pronounced when the 
considerable portion of the wind is accelerated at a distance less than 100 r g 
from BH. Note that the terminal velocity changes form Av 00 /v°° = 0.36 for 
x c = 50 to Av°° /v 00 = 0.19 for x c = 500. The obtained results show that 
the GEF flow can be sufficiently more fast than the flow which is described 
by CAK theory. 

It is illustrative to demonstrate that CAK solution may be obtained from 
( fT7j) by "switching off' smoothly the gravitational redshifting. To obtain the 
continuous transition from the GEF solution to the CAK solution of equation 
fll?]) the following procedure was adopted. We assumed that the radiation 

force in (ITTD is /x ~ dv Idr H — s GM/r 2 The introduced parameter s contin- 

c 

uously changes from 1 (GEF case) to (CAK wind). Numerically calculating 
solutions of (|i~7D for different values of s, applying the inner density bound- 
ary conditions it is possible to demonstrate the continuous transition of these 
solutions from the limiting cases of CAK and GEF solutions. The results are 
shown on Figure 2. 

The introduction of the Paczynski - Wiita ( PW ) ([16]) potential, allows 
to simulate the effects of general relativity. Equations analogous to (|H?D, 
(|25|), (p6|) are derived in Appendix. The results of the numerical integration 
are shown on Figure 4. Introduction of the modified potential can give a 
gain in v°° that varies from Af°° /v°° ~ 0.1 for x c = 30 to Av^/v 00 ~ 0.03 
for x c = 100. 

Approximate relations for a 1? p c . Introducing the following nondim- 
mentional combination e = — , equation (|25j ) will read: 
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\.a\x\ep\ + (a 2 — l)p c — 4a x e = 0. 



(27) 



To obtain an approximate solution of (p7|) we take into account that e << 1. 
Equation ( p7D has two roots: 



(Pc 



From 



8aiea; 2 
28|) we obtain 



a? ± (1 - a\) 



32a|x|e_ 2 ' 
(a? - l) s 



1 + 



(28) 



Pc,l 



1 



,2 : 



Pc,2 



4aia; 2 e 



(29) 



We are looking for the solution with p c > and with a\ = a\/v cr < 1, 
thus the second root of (1297) should be considered: 



p - i — r> ( 3 °) 

where it was supposed that a x << 1. The accuracy of fl3"U|) is Ap c /p c ~ 10~ 3 
(for x c = l,p c = 12.0027, rf ppox = 12.0169). 

To obtain an approximate relation for a\ we substitute p Cj2 from (BS|) to 
(|26|). Taking into account that /3 = Oi/e, the resulting equation reads: 



9 9 a — 1 (1 — a?) 2 a?(l — T) eai . . . 

- 2a x c e 2 + 1 lJ e + lV ; - ^ (1 - a 2 = 0. 31 

a 4ai 2 2a 

Neglecting in (|3l|) terms which contain a^ and a\ e 2 will result in the following 
equation for determining a\. 

2a\(l - V) - 2a\e + = 0. (32) 

a 

The second term in (^) is small compared with other two. It allows to 
obtain the following simple relation for a\. 



Calculating numerically oi, p c , from (|25|), ( f26|) we have found that the value 
of a± depends rather weakly on x c : x c = 50, a\ = 0.034347, p c = 48.0109; 
x c = 500, ai = 0.034348, p c = 0.480112; x c =5000, ai = 0.034369, p c = 
0.004806; (obtained for a = 0.6, T = 0.5). Relation (|33D has a relatively 
good accuracy: for e = 6 • 10~ 5 (B3I) gives al Ppr ~ 0.03430 with accuracy 
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Aai/ai = 10~ 3 . Velocity of the wind at the critical point is approximately 
30 a s (v c = ai l a s ). 

Above the critical point the influence of the d<p/dr (GEF) on the addi- 
tional acceleration of the radiationally driven wind is small compared to the 
effect due to dv/dr. 

3 Discussion 

The most successive model describing winds from hot stars presented so 
far is that of Castor (1970), Castor, Abbott & Klein (1975). Geometry and 
ionization balance differ O-star wind from AGN outflow. Line driven wind 
from O-type star is spherically-symmetric. Outflows in AGN are assumed to 
be originated from the luminous accretion discs and approximately axially- 
symmetric. Powerful X-ray and UV radiation from the disc pose a problem 
of overionization of the outflowing plasma. In Sobolev approximation, which 
is in the background of the CAK theory, the radially streaming radiation 
flux is absorbed in a line transition in a wind with the gradually increasing 
velocity. A photon, emitted by a disc will be red-shifted due to Doppler 
effect. The resultant Sobolev optical depth t x ~ (dv/dr)' 1 . As it was shown 
by CAK, the radiation force that is due to an ensemble of optically thin 
and optically thick lines is proportional to (dv/dr) a . The more dv/dr the 
more effectively a line is shifted to the powerful continuum. In a number 
of studies, CAK theory was applied to AGN in order to explain fast (up to 
~ 0.1c) outflows Arav & Li (1994), Arav, Li & Begelman (1994), Murray 
at al. (1995), Proga et al. (1998), Proga et al. (2000). The CAK theory 
was enhanced by adopting axial symmetry. Ionization balance was studied 
in details simultaneously with 2D hydrodymanical calculations. 

In this paper we developed a theory of winds that takes into account 
effects of the strong gravitational fields. We point out that if a wind is accel- 
erated near super-massive BH a gravitational change of photon's frequency 
must be taken into account. In a strong gravitational field a photon, emit- 
ted by a disc, will be red-shifted due to both the Doppler effect ~ v/c and 
gravitational redshifting Au/u = A0/c 2 . We argue that taking into account 
gravitational redshifting can substantially change the wind dynamics and 
structure. Although it should be mentioned that the developed theory (in 
the adopted toy model) cannot be directly applied to explain outflows from 
close to BH. An exact self-consistent solution of this problem is possible 
only in general relativity. We avoid this sophisticated task by considering all 
equations in flat space and time. 

Basing on considerations of Sobolev we conclude that the greater d(f)/dr 
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the more effectively a line is shifted to the extent where the radiation flux 
is unattenuated. In such a case gravitational field "exposes" a wind to the 
unattenuated continuum. Note that this effect is independent of the wind 
dynamics, it works also in a medium with v — 0. We determine a wind 
accelerated in this regime: " Gravitationally Exposed Flow" (GEF). 
The main goal of this paper has been to compare a solution for GEF with 
that obtained from a standard line driven wind theory ( CAK ). To solve 
this problem a very simple input physics was assumed: spherical symmetry, 
constant temperature and no ionization balance. We found that in such a 
case the radiation force /; ~ [dv/dr + ^d(f> / 'dr)~ l . Numerical analysis shows 
that the introduced effect changes drastically a slope of a solution curve at 
the bottom of the wind. It is clearly a result of the fact that the role of the 
gravitational field is important when the field is high an velocity gradient is 
low. As soon as dv/dr becomes sufficient, Sobolev effect efficiently stimulates 
the acceleration of the wind. Obtained equation of motion is nonlinear with 
respect to dv/dr but this nonlinearity is different form that of CAK. Critical 
point of the equation of motion is neither a sonic point nor a CAK critical 
point. The position of the critical point r c is the only free parameter of 
the problem. Adjusting r c in order to satisfy density boundary conditions 
we compared solution for GEF with that of CAK which was obtained for 
the same physical and boundary conditions. Numerically solving equations 
describing spherical symmetric, stationary outflowing wind we found that 
gravitational redshifting can make acceleration up to 35% (for r cr = 50 r g ) 
more efficient. 

In order to take approximately into account effects of general relativity, 
we made use of modified potential (Paczynski-Wiita potential). Numerical 
analysis demonstrated that if the critical point is located as far as 30 r g 
100 r g the difference between PW case and Newtonian case is Av^/v 00 = 
0.1 ^0.03 

Acknowledgements Author is deeply grateful to G.S. Bisnovatyi-Kogan 
for continuous and stimulating attention to this work. 

4 Conclusions 

1. A theory of winds accelerated by the radiation pressure in lines with 
account of gravitational redshifting of photons is developed. A system 
of equations describing stationary, spherically-symmetric, isothermal 
flow is derived. 

2. A solution of these equations is obtained numerically for two cases: 
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a standard line-driven wind (CAK theory), Gravitationally Exposed 
Flow (GEF)- a wind that is accelerated by the radiation pressure in 
lines if to take into account gravitational redshifting. It is shown that 
an increase of up to 35% in v 00 can be obtained. 

3. To take approximately into account effects of general relativity, Paczyn- 
ski - Wiita potential is adopted. A wind solution is calculated for this 
type of potential and comparative analysis with Newtonian case is pre- 
sented. 

4. The developed theory can be used to explain fast outflows from AGN. 
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5 Appendix 



Substituting Paczynski - Wiita potential into equation (|T|) will obtain equa- 
tion of motion: 



a \ dv 2a 2 GM 



La 



dr r (r — r g ) 2 47rr 2 c 



Lak 



Air 



4%r 2 c \ \M\av th , 



'dv 1 GM 



vr - — | — 



dr c (r — r g y 



(34) 



Scaling all variables according to (|T8|), omiting tilda, equation (Rl) reads: 



F(p, v, x) 



x 



P 



V I 

P 2 \ P + 2(x- l) 2 



0. 



(35) 



where p = dv/dx, /i is determined by (|2l|) and a±, (3 by ([201). Making use 
of (^) and ( p3|) will result in the following relations: 



a 1 /-I 2\ 2 o 2 i 

-(1— a x )x c p c — 2a 1 x c + — 



xt. 



a 



2 \x c -lf 



-V)-{l-a{) 



2^ P 



x,. 



2a (x c - l) 2 



c 2 p c + cip c + c = 0. 
Coefficients determined from the following relations: 



0, (36) 
(37) 



C-2 

Co 



2a 1 x c , 



\)x\P 



2(x c - l) 2 ' 

-x c (P - 1)0 + a 2 {-2{x c - If - f3x c 



Equation for fi reads: 



a 



- -. z. 

J 2 



ix c - 1)* 



Pc + 



p 



2(Xr.-l) 



1-a 



(38) 
(39) 

(40) 



(41) 



Equation (41) is analoguis to fl24]) except the last term in brackets. For a 
given x c , equations (|36|), fl37j ) are solved numerically to determine p c and a\. 
Then ([35]) is integrated numerically as described in the text. Comparative 
results of the numerical solution of fl3"5]) are depicted of Figure 4. 
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Figure 1: Solution of the wind equation (19|). Solid line - GEF solution, 
dashed line - CAK solution. Crosses indicate GEF critical points, circles - 
CAK critical points. Curves for x = 8000 are graphically indistinguishable. 
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Figure 3: Density distribution in the wind. The solution curves calculated 
for x<Z EF = 50, x G c EF = 75. 
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Figure 4: Solid line indicates GEF solution for the Newtonian potential, 
dashed line GEF solution for Paczynski - Wiita potential. 
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